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ABSTRACT 
We present the discovery of a super-earth-sized planet in or near the habitable zone of a sun-like 
star. The host is Kepler-69, a 13.7 mag G4V-type star. We detect two periodic sets of transit signals in 
the three-year flux time series of Kepler-69, obtained with the Kepler spacecraft. Using the very high 
precision Kepler photometry, and follow-up observations, our confidence that these signals represent 
planetary transits is >99.1%. The inner planet, Kepler-69b, has a radius of 2.24lg'29 -R©aiid orbits 
the host star every 13.7 days. The outer planet, Kepler-69c, is a super-Earth-size object with a radius 
of l-7lo'23 i?®and an orbital period of 242.5 days. Assuming an Earth-like Bond albedo, Kepler-69c 
has an equilibrium temperature of 299 ± 19 K, which places the planet close to the habitable zone 
around the host star. This is the smallest planet found by Kepler to be orbiting in or near habitable 
zone of a Sun-like star and represents an important step on the path to finding the first true Earth 
analog. 

Subject headings: planetary systems; stars: fundamental parameters; stars: individual (Kepler-69, 
KIC 8692861, KOI-172); techniques: photometric; stars: statistics 



1. INTRODUCTION 

NASA's Kepler mission was launched in 2009 with the 
primary aim of determi ning the abundance of Earth-size 
planets in our Galaxy (jBorucki et al.lf2010( ). The 0.95- 
m optical telescope is in an Earth-trailing heliocentric 
orbit, providing a stable thermal environment for the 
production of very high precision photometry. It com- 
prises 42 CCDs wi th a combined field of view of 115 dcg^ 
ijKoch et al.ll2010f ). The brightness of over 160,000 stars 
is measured in 29.4-min integrations, and these data arc 
searched for transiting planets, which appear as periodic 
dimming in brightness. 

While the first true Earth-analog, a planet the size of 
the Earth and in the habitable zone of a Sun-like star, 
has yet to be discovered, several important milestones 
towards achieving this goal have been met using data 
from the Kepler spacecraft. The first planet found by 
Kepler to be orbiting in t he habitable zone o f a Sun- 
like star was Kepler-22b (jBorucki et al.ll2Q12D . With 
a radius of 2.4 Rq, Kepler-22b likely does not have a 
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rocky surface and is po ssibly a temperate ocean world 
([Rogers fc Seaged 12010( 1. Kepler h as also been used to 
discover planets th at are Earth-siz e (iBatalha et al.l 



Fressin et a l. '2012') and smaller (jMuirhead et al.l 



2011 



2012 



Barclav et al. 2013) but these planets are likely to be too 



hot to host liquid water at their surfaces. The three cat- 
alogs of planet candidates pu blished by the Keple r mis- 
sion (jBorucki et all l2011allbl : IBatalha et all 120131 ) have 
detailed over 2300 candidates, and have each contained 
progressively smaller candidates in longer orbits. How- 
ever, the observation timespan searched in each catalog 
(43 days, 137 days and 1.5 years, respectively) have thus 
far been too short to allow for three transits (necessary 
for period confirmation) of a potentially rocky planet in 
the habitable zone of a sun-like star. 

Kcplcr-69 has been observed for nearly three years by 
the Kepler spacecraft (observational quarters Q1-Q12) 
in both long cadence mode (continuous 29.4-min inte- 
grations), and during Q3-Q6 and one month of Q7 in 
short cadence mode (continuous 58.8-sec integrations). 
Transit-like signals in the flux time series, occurring ev- 
ery 13.7-days wi th a depth of 600-ppm, have been previ- 
ously reported (JBorucki et al.ll2011bl : iTenenbaum et al.l 
l2012t IBatalha et al.ll2013() . This signal was classified as 
a Kepler Object of Interest (KOI) with catalog number 
KOI-172. 01. No other planet candidates were detected 
in the first 1.5 years of observations. Here we report 
the detection of a second planet candidate, KOT172.02, 
with a period of 242.5 days, which was detected in a 
search of the first 2 years of observations. Since plan- 
ets of this size and period cannot currently have their 
masses measured (by, for instance, radial velocity mea- 
surements) and their planetary nature thus confirmed, 
we rely on statistical analysis of the likelihood that KOI- 
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172.01 and KOI-172.02 are planets for validation of their 
planetary status. 

2. STELLAR PROPERTIES 

Kepler-69 appea rs in the Kepler Input Catalog (KIC, 
iBrown et aLll2011[) as KIC 8692861 and has a magnitude 
of 13.7 in the Kepler bandpass. To derive stellar proper- 
ties we obtained a high resolution spectrum of Kepler-69 
using the HIRES spectrograph on the Keck I telescope 
on 23 June 20 11 using the setu p of the California Planet 
Search group (jMarcvet al.l [20Qi). 

The stellar properties of Kepler-69 were determined by 
performing a x^ fit of the spectrum with a library of 750 
observed spectra of F-M type stars that have accurate 
parallaxes. The library spectra have effective tempera- 
tures spanning, Teg = 3500K-7500 K and surface grav- 
ities, log 5 = 2.0-5.0. The weighted mean of the ten li- 
brary spectra with the lowest x^ values were adopted as 
effective temperature, stellar surface gravity and nietal- 
licity of Kepler-69. The weights used were inversely pro- 
portional to the Pythagorean distance of stellar parame- 
ters from the median value, leading to an iteration that 
converges quickly. Choosing 10 library matches from 
which to determine the weighted average offers a statisti- 
cal buffer against outliers. Best-fitting stellar properties 
are provided in Table [TJ ^^^^^_^^^^^_^_^ 

We also ran an SME analysis (jValenti fc Fischeil[2005[) 
of the Keck spectrum of Kepler-69. The SME results 
were in excellent agreement our spectrum matching tech- 
nique, albeit with l ower uncertainties. However, fol- 
lowing [Torres et al.l ()2012D we inflated the uncertaities 
derived from the SME analysis to account for system- 
atic biases. These inflated uncertainties were similar 
to those found from the spectral matching method de- 
scribed above so we chose to use spectral matching de- 
rived parameters in the remaining analysis. 

We searched Kepler short cade nce data for solar- 
like oscillations using the method of iHuber et al.l ()2009l ) 
but did not detect any. S ince oscillation am plitudes 
scale with stellar luminosity (jKieldsen fc Bedd ing 1995), 
evolved subgiant and giant stars show oscillations that 
are readily detectable with Kepler data of 13.7 mag 
stars. The non-detection of oscillations confirms that 
Kepler-69 is a dwarf star with a lower limit of log g > 3.8, 
consistent with the stellar properties derived from spec- 
troscopy. 

We matched the spectroscopically derived temper- 
ature, surface gravity and metallicity to a fine grid 
of evolutionary model s from the BASTI database 
(|Pietrinferni et ani2004D in order to estimate the star's 
mass and radius. We report the median and the 
1-cr region of each parameter derived from Monte- 
Carlo simulations using the spectroscopic values in Ta- 
ble[Tl Consistent pararn eters are found usin g Yonsei-Yale 
(iDemargue et all 12004! ) and Dartmouth ([Potter et al.l 
|2008() evolution models. Kepler-69 is a main-sequence 
G4V star, with a somewhat lower mass and metallicity 
than the Sun. 



TABLE 1 
Stellar properties of Kepler-69 



3. TRANSIT ANALYSIS 



3.1. 



Transit detection by the Kepler pipeline 

The Transiting Planet Search (TPS, iJenkins et al.l 
2010cl: iTen enbau met all 1201 2[ ) and Data Validation 
(DV, iWu et al.. ,20l3) modules of the Kepler Pipehne 



Parameter 



Value 



Teff (K) 

logg (dex, cgs units) 

Metalicity [Fe/H] 

Mass (Mq) 

Radius (Rq) 

Luminosity (Lq) 

Density (g cm~^) 

Radial velocity zcropoint (km s^-*^) 



5638 ± 168 
4.40 ±0.15 
-0.29 ±0.15 

o.8iotg:»? 

0.93l«;i« 

-1+0.37 



0.80] 
1.37^ 



0.22 
0.81 



-0.55 
-38.7 ±0.1 



(jJenkins et al.l l2010bD identified two transit-like signa- 
tures in a search of Ql-QlO data. The first signature 
was due to KOI-172.01, a planet c andid ate previously 
report ed by iBorucki et al.l (j2011b[) and IBatalha et al.l 
(|2013() . The second signature was new and had a pe- 
riod of 242.4579 ± 0.0056 days and a depth of 258 ± 19 
ppm and has been designated KOT172.02 (Burke et al. 
in prep.). The planet radius provided by the Q1-Q12 
DV report for KOI-172.02 was 1.4 i?®. A signal-to-noise 
ratio (SNR) of 73 was calculated for the transits of KOI- 
172.01 and 15 for KOI-172.02. The SNRs for both planet 
can didates are si gnificantly above the formal threshold of 
7.1 (jJenkins et a l. 2002,), giving us very high confidence 
that neither of these detections are due to random or 
correlated noise. 

3.2. Light curve preparation 

There are instrumental features near two of the KOI- 
172.02 transits. One was a su dden pixel sensitivit y 
dropout due to a cosmic ray hit (IStumpe et al.ll2012aD, 
the o ther is a coronal mass ejection (j Thompson et al.1 
|2012[ ). Both apertu re photometry and pr e -search data 
cond itioned (PDC, IStumpe et all 1201 2al: iSmith et all 
120121 ) flux time series data stored at the Mikulski Archive 
for Space Telescopes (MAST) insufficiently corrected the 
out-of-transit flux levels around the instrumental signals. 
This led to a measured transit depth that was under- 
estimated by 20%. We instead chose to use a devel- 
opmental version of the PDC error-corrected data that 
relies on error correction via wavelet-based band split- 
ting (jStumpe et al.|[2012bD . This newly corrected data 
preserves stellar variability. We removed stellar variabil- 
ity using a nu merically efficient discreet cosine transform 
(|Garciall2010l ) to implement a non-parametric penalized 
least squares smoothing of the fiux time series. Data 
from individual quarters were normalized and then com- 
bined. Transits were treated as missing when smoothing 
the data. 

3.3. Transit fitting 

We performed a fit of both pla n et can didates simul- 
taneously using a iMandel fc Agoll (|2002[ ) transit model 
with quadratic limb darkening. Limb darkening pa- 
rameters were c omputed by interpolating the tables of 
IClaret fc Bloemcn (2011) and kept fixed. 

The parameters included in our fit were: mean stellar 
density, photometric zeropoint, and for each planet the 
transit epoch, orbital period, the planet to star radius 
ratio, impact parameter, and eccentricity vectors esinw 
and ecosw, where e and w are the eccentricity and ar- 
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Fig. 1. — The Kepler light curve of KOI- 172. The upper panel 
shows the transits of KOI-172.01 folded on the orbital period and 
the lower planet shows the transits of KOI-172.02. The blaek points 
show the observed data and the blue points are binned data. In the 
upper panel 100 observed points are included in each bin and in 
the lower panel 12 point are included in each bin. The uncertainty 
on the binned data is the standard deviation divided by the square 
root of the number of point included in the bin. The best fitting 
transit models are shown overplotted in red. 

gument of periastron of the orbi t, respectively. We uti- 
lized the emcee implementation (|Foreman-Mackev et al.l 
|2012[ ) of an affinc-invariant Markov C hain Monte Carlo 
(MCMC) algorithmf Good man fc We arc 2010) to calcu- 
late posterior distributions for these parameters. The 
mean stellar density measured from asteroseismology was 
used as a Gaussian prior, all other parameters were given 
uniform priors except for the two eccentricity vectors. 
Using the e sin w and e cos uj parameterization results in 
an implicit linear prio r in e (|Fordll2006HBurke et al]l2008l : 
lEastman et al.|[2013D . We correct for this by enforcing 
a 1/e prior on eccentricity. Our likelihood function is 
therefore 



C ~ exp 




[Pa - PmY 



Ine 



(1) 



where pa is the mean stellar density derived from astero- 
seismology and Pa is the model mean stellar density, x^ 
is the usual sum of squared deviation from the model 
weighted by the variance. 

The median and central 68% of the posterior distribu- 
tion (equivalent to the 1-cr uncertainty) of each parame- 
ter arc shown in Table[2] Uncertainty in the stellar radius 
was included as a component of these calculations. This 
analysis yielded planet radii of 2.24^Q2g and 1.7i;to'23 
i?e for KOI-172.01 and KOI-172.02, respectively. 



TESTS TO ASSESS THE VALIDITY OF INTERPRETING 
KOI-172.01 AND KOI-172.02 AS PLANETS 

4.1. Tests on the Kepler flux times series data 



The Kepler photometric data for KOI-172.01 and KOI- 
172.02 were examined for evidenc e suggestive of a fals e 
positive in the manner described in lBatalha et al.l()2013D . 
No such evidence was found and hence both candidates 
were afforded a place within the most recent Kepler 
planet candidate lislO (Burke et al. in prep.). Specif- 
ically, odd and even numbered transits have consistent 
depths; there is no sign of a secondary transit; and there 
is no significant shift in the phot o-centroid of the star in 
transit relative to out of transit (iJenkins et al.|[2010aD . 

The Data Validation (DV) module of the Kepler 
pipeline creates three images; an average out-of-transit 
image taken from nearby but not during the transit 
events, an average in-transit image and a difference im- 
age which is the difference between the out-of-transit im- 
age and the in-transit image. A fit of the pixel response 
function (PRF) to the difference image gives the posi- 
tion of the transit source while a fit of the PRF to the 
out-of-transit image yields the position of the target star. 
The offset between the difference image position and the 
target star position can be used to identify false positive 
scenarios. Conversely, if no significant offset is detected, 
the uncertainty in the difference image centroid position 
can be used to calculate a radius of confusion outside of 
which a false positive source can be excluded. No signif- 
icant offsets are measured between the difference image 
and target position for KOI-172.01. The 3-cr radius of 
confusion was measured by DV at 0.16 arcsec. Only two 
transits were included in the PRF centroid offset metric 
in the DV report for KOI-172.02 owing to data artifacts 
and transits of KOI-172.01 occurring near the transits of 
KOI-172.02. We wished to use all five transits of KOI- 
172.02 in our centroid analysis so we created difference 
images in the same manner used by the Kepler pipeline 
(jBrvson et al.l [201 30 but included all observed transits 



of KOI-172.02. The significance of the offset from the 
remaining images was less than S-g which is used as 
the w arning threshold (jBatalha et al.|[20T3l : iBrvson et al.l 
I2013D . The radius of confusion was estimated to be 1.5 
arcsec. Given the above described tests, there is no rea- 
son to suspect KOI-172.01 or KOI-172.02 are false posi- 
tives based on the Kepler data. 

We measured a transit depth for KOI-172.02 of 350 
ppm from the MCMC analysis. The faintest a star could 
be and still produce a transit of this depth can be calcu- 
lated if we a ssume a total eclipse of that background star 
(jChaplin et al. 2013|). Under this assumption the faintest 
a false positive star could be is 7.9 mag fainter than 
Kepler-69. This is the limit we use in all false positive 
calculations for KOI-172.02. For KOI-172.01 the transit 
depth is 597 ppm, which sets the maximum brightness 
difference between the target star and a false positive of 
7.3 mag. 

4.2. Imaging of Kepler-69 

We obtained optical F-band images of Kepler-69 us- 
ing the Nickel 1-m telescope at Lick Observatory and 
searched for nearby stars. No stars are seen 2-5 arc- 
sec from Kepler-69 brighter than 19th magnitude; stars 
closer than 2 arcsec could not be ruled out. Kepler-69 

^ The Kepler planet catalog and DV re- 
ports are hosted at the NASA Ex oplanct Archive, 
|http : //exoplanetar chive ■ Ipac . caltech . edu/ | 
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TABLE 2 








Parameters 


FROM MCMC ANALYSIS 






Parameter 






Kcpler-69b {KOL172.01) 


Kepler-69c (KOL172.02) 


Stellar density (g cc^"^) 






1.05 


fU.4B 
-0.27 




-RpZ-R* 






0m207ttSS,Ts 




0.0168t°:Sg 


Mid-point of first transit 


(BJD- 


- 2454833;!^ 


137.8414«;««i« 




150.870lHi^ 


Period (days) 






1Q ^000.11+0.000035 




242A613t'ofoll 


Impact parameter 






o.i5lg:iS 




n 19+0-21 


ecostj 






omt°oil 




-0.01+°!^ 


esinw 






o 07 + 0.09 




-0.02l°;?« 


Planet radius (i?®) 






r, r,<+0.44 




1-711°:^^ 


a//?* 






21.8+^-? 




148l?o 


Scmimajor-axis (au) 






0.094lH?^ 




0.64lg:l? 


Inclination (deg) 






89.62+E;-^® 




89-85t°:°| 


Eccentricity 






n ifi+oi'^ 




o.i4tg:ig 


Equilibrium Temperature 


'.{Kt 




7791- 




299t^g 


Limb darkening parameters 




{0.4012, 


0.2627} 



''BJD-2454833 is the standard Kepler time system, known informally as BKJD. 

''Equilibrium temperature assumes an albedo of 0.3 and the thermal recirculation parameter of 1.0. 



has a y-band magnitude of 14.0 (lEverett et al.ll2012l ). 
therefore we cannot exclude stars between V=19-21.9 
for KOI-172.G2 and V=19-21.3 for KOI-172.01 based on 
the Lick image and the transit depths. 

There is a UK Infrared Telescope (UKIRT) J-band im- 
age of a field containing Kepler-69. The closest star is 
3.2 arcsec west of Kepler-69 and is 9 mag fainter. This 
close-by star can be ruled out as the source of the transit 
signal from either planet candidate based on the PRF 
centroids. A magnitude dependent radius of confusion 
around Kepler-69 was calculated for use in false positive 
analysis. This confusion radius was converted to the Ke- 
pler m agnitude system (jBrown et al.ll20lit iHowell et al.l 
|2012| ) and is shown in blue in Figure [H 

4.3. Looking for a second star in HIRES spectra 

The slit width used to collect the HIRES spectrum 
of Kepler-69 was 0.87 arcsec. This allowed us to search 
the HIRES spectrum for additional stars within 0.43 arc- 
sec of the target. We calculated the best-fitting two- 
spectrum model using observed spectra with T^s rang- 
ing from 3500-6000 K. We searched over a range of sec- 
ondary star radial velocities and flux ratios relative to 
the target star. No second star was seen in the spectra 
down to the confusion limit of AKp = 4.0. We then 
injected a range of observed spectra into the spectrum 
of Kepler-69. Given we measure a radial velocity for 
Kepler-69 of —38.7 ±0.1 km s~^, virtually all sufficiently 
bright background stars would have been detected. G 
and K dwarf physical companions to Kcpler-69 that are 
further than approximately 20 au could evade detection 
owing to them having a radial velocity that we could not 
distinguish from Kepler-69. Companions with spectral 
type later than MO are detected at any distance from 
the target (subject to falling into the spectrograph's slit) 
provided they are brighter than AKp = 4.0. 

5. ON THE PROBABILITY THAT KOL172.01 AND 
KOL172.02 ARE FALSE POSITIVES 

In this section we first consider the probability that 
KOI-172.01 and KOI-172.02 are in fact transits or 
eclipses of a background star. We find that the chance 



of this occurring for KOI-172.01 is 0.01% and for KOI- 
172.02 it is 0.04%. We are able to rule out the transit- 
like signals being due to a heirarchical triple stellar sys- 
tem and therefore validate KOI-172.01 and KOI-172.02 
as confirmed planets. 

We then go on to look at whether either of the two 
planets orbit a physical companion star. While for KOI- 
172.01 the probability of this is 0.01%, for KOT172.02 it 
is <1% and hence cannot be definitively ruled out, but 
is unlikely. 

5.1. False positive scenarios considered 

Given that the KOI-172.01 an d KOI-172.02 detect ions 
are unlikely to be due to noise () Jenkins et al.|[200^ . we 
can safely assume that they either represent planets or 
astrophysical false positives. Scenarios that could cause a 
transit-like signal that are not caused by a planet orbiting 
Kepler-69 are: 

1. Stellar binaries with grazing eclipses (the impact 
parameter >1). 

2. Background eclipsing binariefl 

3. Background star-planet systems. 

4. An eclipsing binary physically associated with 
Kepler-69 in a heirarchical stellar triple system. 

5. A planet orbiting a stellar physical companion of 
Kepler-69. 

We are immediately able to dismiss scenarios 1 and 4 
because the best-fitting models representing these sce- 
narios fit the observed data worse than the best-fitting 
transiting planet model at a level >3cr. This is because a 
stellar eclipse across Kepler-69 or a physical companion 
would be more V-shaped than we observe. Additionally, 
these scenarios are predicted to be very uncommon for 

^ The phrase background stars is used to refer to stars that 
appear fainter than Kepler-69 and are not physically associated. 
They may be cither in the foreground or background with respect 
to Kcpler-69. 
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Fig. 2. — The regions of parameter spaee where a false positive source can be excluded are shown colored. The exclusion region from 
the transit depth is shown in red, from the UKIRT imaging in blue and from the centroids in cyan. The white hashed region shows the 
constraints from the Keck spectral analysis. In this region contaminating stars are only allowed if they have a radial velocity within 10km 
s~^. The region in white is not excluded and must therefore be accounted for in our false positive analysis. 



sub- Jupiter-sized planet candidates (jFressin et al.ll2013h . 
The remaining false positive discussion is restricted to 
scenarios 2, 3 and 5. 

5.2. Background false positive scenarios 

To assess the probability that the transit-like signals 
we detect arc not associated with the Kcplcr-69 system 
but are on a background star, we simulate the stellar 
background population that could be responsible. We 
simulated the number of stars within 1 deg of Kepler-69 
using the Bcsancon Galaxy model with kinematics en- 
abled (.Rob in ct al. 200^. The Galactic model predicts 
there are 432,118 stars, of which 44,185 fall in the bright- 
ness range 14.0 < V < 22.0 and are not excluded by the 
secondary star test performed on the Keck spectrum. We 
integrated the area shown in white in Figure [5] to calcu- 
late the undetected stellar population in the background 
of KOI-172.02. The estimated number of stars in the 
background is 0.018. If we now consider that only 2.6% 
of stars observed by Kepler host either a transiting planet 
candidat e or are an eclipsing binary, excluding c ontact 
binaries (jBatalha et al.ll2013nSlawson et al.ll201lL Burke 
et al. in prep.), the predicted number of background 
stars which could cause a false positive is 0.00046. 

We used a Bayesian appr oach similar to that applied 
by iFressin et ahl (j2012[) and jBarclav et all (|2013[) to the 
validation of the planets orbiting Kepler-20 and Kepler- 
37 to quantify the false positive probability. This tech- 
nique compares the a priori chance of finding a planet the 
size of KOI-172.02 (the planet prior) to that of finding a 
background false positive. The ratio of the planet prior 
to the false positive probability is used to calculate our 
confidence in the planet interpretation. The population 
of super-earth size planets orbiting with periods longer 
than 85 days is not well constrained. We assume that the 
super-earth occurrence rate is the same at long orbital 



periods as it is at shorter orbital periods when period is 
measured in logarithmic units. There is evidence that 
the occurrence rate of super-Earth-size planets is even 
high at longer orbital periods than at short periods. For 
example iGassan et al.l ()2012f ) find 62^37% of stars host a 
su per-Earth-mass plan et. Therefore, using the statistics 
of lFressin et al.l (|2013f ) who derived a super-Earth occur- 
rence rate of 23±2.4% may well lead to an underestimate 
of our planet prior. 

There are 666 super-e a.rth-size (1.25- 2 .0 Rm ) planet 
cand idates tabulated in iBatalha et al.l ()2013() . How- 
ever, IFressin et al.l (|2013D predict that 8.8% of these are 
false positives which provides an estimated number of 
super-earths found by Kepler of 607. The catalog of 
IBatalha et al.l (J2013D was based on continuous (Q1-Q6) 
observations of 138,253 stars. Therefore, the occurrence 
rate of transiting super-Earth-size planets is 0.0044 per 
star. We can compare this to the background/foreground 
false positive population which is 0.00046, as explained 
above. If KOI-172.02 were the only planet candidate as- 
sociated with Kepler-69, the likelihood of it being a true 
planet would be 0.0044/(0.00444-0.00046) = 0.91 . How- 
ever, KOI -172.02 is in a two pla net candidate system. As 
shown bv lLissauer et al.] (|2012D . having multiple transit- 
ing planet candidates associated with a star increases 
the probability that these candidates are real planets by 
a factor of '^15 compared to single planet candidate sys- 
tems. This multiplicity boost increases our confidence 
that KOI-172.02 is a bona fide planet physically associ- 
ated with the Kepler-69 system to 99.3%. 

We note that Kepler-69 has a lower metallicity than the 
Sun and for large planets a low metal cont ent has been 
shown to hinder planet occurrence (e.g. IMavor et al.l 
l2011t ). However, for there appears to be no link between 
plane t occurrence and stellar metalicity (jBuchhave et al.l 
I2OI2I ) for transiting super-Earth-size planets. Therefore 
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we do not factor the stellar metalicity into our our false 
positive calculations. 

The centroid confusion radius for KOI- 172. 01 is much 
smaller than for KOI-172.02 (as shown in Figure [2|). Us- 
ing the approach outlined above, our confidence is 99.9% 
that KOI-172.01 is a planetary mass body physically as- 
sociated with the Kepler-69 system. We can therefore 
conclude that the signals we observe are due to bona fide 
planets. The inner and outer planets are named Kepler- 
69b and Kepler-69c, respectively. 

5.3. The possibility that the planets in the Kepler-69 
system orbit a physical companion. 

While the above analysis leads us to conclude that 
Kepler-69b and Kepler-69c are both sub-stellar bodies 
physically associated with the target system, this does 
not exclude a scenario where one or both of the planets 
orbit a fainter stellar companion to Kepler-69. 

We performed a series of MCMC transit analyses with 
no prior on the stellar density with the aim of constrain- 
ing the range of stellar types the planets could orbit. 
First we assumed that both planets orbit the same star. 
If this is the case then this star would have a stellar den- 
sity of 1.04to;3i g cm~'^ with a 3-tT upper limit of 3.6 g 
cm~^. This rules out both planets orbiting a star with a 
radius less than 0.64 Rq, wh ere the radius is ca lculated 
using Dartmouth isochrones (jDotter et al.ll20d8l ) and as- 
sumes Kepler-69 and stellar companion co-evolved. We 
tried including a fixed dilution in the transit depth of 
98% to simulate both planets being around a star l/50th 
the brightness of the primary star in the system. In this 
case we found a 3-a upper limit on the stellar density 
of 2.8 g cm~'^. However, a companion hosting a transit 
that who's light is diluted by 98% should have a density 
of 9.7 g cm~'^ based on isochrones, therefore we can rule 
out a scenario where both planets orbits a much fainter 
companion. If both planets orbit a companion, it must 
be a G or early K-type star. 

We now consider the case that both planets orbit dif- 
ferent stars. Fitting only the transits of Kepler-69b gave 
a stellar density upper limit of 6.2 g cm~'^ which would 
place Kepler-69b around a star with a radius of at least 
0.48 Rq. This rules out Kepler-69c orbiting a star with a 
spectral type later than K9V. Fitting just for Kepler-69c 
implies a star with a mean stellar density of l.Stio § 
cm~^ and the 3-a upper limit on density is 8.7 g cm~^. 
Including a dilution in the transit depth of 98% loosened 
our 3 — CT upper limit on the stellar density to 15.5 g 
cm~^. A star with a density of 15.5 g cm^'^ that co- 
evolved with Kepler-69 would have a radius of 0.31 Rq 
- an M2-type star. Therefore, in the case of a physical 
companion, Kepler-69c cannot orbit a star cooler than 
an M2-type. 

To quantify the probability that either of the plan- 
ets orbit a companion star, we constructed three Monte 
Carlo simulations to estimate the properties of a stellar 
companion with different assumptions: both planets or- 
bit a companion to Kepler-69; Kepler-69b orbits Kepler- 
69 and Kepler-69c orbits a companion; and Kepler-69c 
orbits Kepler-69 and Kepler-69b orbits a companion. 

In the simulations we assume the orbital period is de- 
scribed by a log- normal distribution with mean logP = 
5.03 and standard deviation CTiogp — 2.28 where period 



is in days (jRaghavan et al.|[2010D and eccentricity is uni- 
form from zero to one. The primary to secondary mass 
ratio q is assumed to be uniformly dist ributed (this ap- 
proxi mates the distribution found by iHalbwachs et all 
|2003| ). The argument of periastron is taken to be uni- 
formly distributed. We ran the simulation 10^ times. 
At each loop in the simulation we calculated the radius, 
density and bright ness of the companion based on Dart- 
mouth isochrones (jDotter et a l. 2008). We calculated the 
binary separation from the orbital period and mass ratio. 
The radial velocity at the point of observation of the sec- 
ondary was determined by calculating a radial velocity 
curve and selecting the time of observation at random. 

We then calculated the proportion of secondary scenar- 
ios from our simulation that are unphysical or excluded 
based on observations. Unphysical scenarios are those 
where the binary separation is to o small to allow a sta- 
ble planetary orbit of 24 2.4 days ()Rabl fc Dvoraklll988l: 
iHolman fc WiegertJ[T999h . We used the diluted MCMC 
transit analyses' 3-ct upper limit on stellar density to ex- 
clude companions that arc too dense or faint. Using the 
constraint from the Keck spectrum, if a companion had 
a brightness within 4 mag of Kepler-69 and a radial ve- 
locity difference between it and the primary of >10 km 
s^^ (the limit from our analysis of the Keck spectrum) 
it was excluded. Companions that would be seen by the 
J-band image or would induce a detectable centroid shift 
were also excluded. We find that a companion star is able 
to host Keplcr-69b less than 0.01% of the time. There- 
fore, we are confident that Kepler-69b orbits the target 
star. However, in 1.9% of the simulations Kepler-69 has 
a stellar companion capable of hosting Kepler-69c. If 
we assume that in scenarios with a viable planet host- 
ing companion, half the time the planet orbits the pri- 
mary and half the secondary, then our confidence that 
Kepler-69c orbits the primary is 99.1%. However, no 
multiplicity boost is included in the above calculation. It 
is unclear whether a companion is more likely to host 
transiting planets if the primary hosts transiting planets 
compared to a field star. However, if we were to apply 
this boost, the probability that Kepler-69c orbits Kepler- 
69 is 99.93%. The true probability that Kepler-69c orbits 
a companion star lies somewhere between the 0.06% and 
0.9%. 

While we cannot conclusively rule out the possibil- 
ity that Kepler-69c orbits a smaller star and therefore 
is larger and cooler, the chances are low. With more 
follow-up observations we may be able to further restrict 
the allowed parameter space and increase our confidence 
that Kepler-69c orbits the primary star. Specifically, con- 
tinued non-detections in future radial velocity observa- 
tions will increase the size of the region around Kepler-69 
where we can rule out any companion star, while a deep 
adaptive optics observations can help rule out compan- 
ions further from the star than around 0.5 arcsec. 

In the remaining analysis we will assume that Kepler- 
69c orbits the target star. 

6. ON THE COMPOSITION AND TEMPERATURE OF 
KEPLER-69C 

It is not clear where the boundary between the ra- 
dius of a rocky (or terrestrial) planet and that of an 
giant volatile rich planet lies. There is very likely not 
an absolute boundary but a range of cutoffs which de- 
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pend on the initial comp osition of the proto-planetary 
disk () Valencia et alj [20071) . In any case, the l-tr uncer- 
tainty we determine for the radius of Kepler-69c makes 
predicting its composition difficult. The l-cr lower bound 
of 1.48 i?0 may well represent a rocky planet but the 
upper bound of 2.04 probably represents a fairly volatile 
rich planet. If Kepler-69c were volatile rich is may be 
a water world and quite unlike any planet in our Solar 
System. 

To estimate the equilibrium temperature of Kepler-69c 
we must make some assumptions about the albedo and 
the thermal energy redistributi on of the plane t . We use 
the same equation as used in iBatalha et al.l (120131 ) to 
calculate equilibrium temperature, 



Tcq-^eff(i?,/2a)l/2[/(l-^B)]'/^ 



(2) 

Teff is the effective temperature of the star, _R^ is the stel- 
lar radius, a is the planet's semi-major axis, / is the ther- 
mal redistribution parameter (where f ~ 1 refers to full 
thermal recircularizatioii], and Ab is the Bond albedo. 
In the Batalha et al. (20131) Kepler planet candidate cat- 
alog they assume f — 1 and As ~ 0.3 (the same albedo 
as Earth). Using these values results in a temperature 
of 299^^9 K, significantly lower than the boiling point of 
water. If we assume a lower Bond albedo of 0.1 (similar 
to that of Mercury) the temperature of the planet in- 
creases to 318J;2i K, still in a regime where water could 
exist in a liquid form (neglecting the effects of an atmo- 
sphere). In Figure [3] we show the range of equilibrium 
temperatures that Keplcr-69c would have as a function of 
Ab- With a reasonably thick atmosphere and a rotation 
rate similar to Earth, Kepler-69c would have a thermal 
redistribution parameter of / ~ 1 (for reference Earth 
has a value of / = 1.1 and for Venus / ~ 1). 

An atmosphere can dramatically increase t he su rface 
temperature of a planet (c.f. Venus). iKasting f|2011|) sug- 
gests the upper boundary on the temperature of a hab- 
itable planet lies around 309 K - Kepler-69c likely falls 
below this limit i f it has a simfiar o r higher albedo than 
Earth. However, iSelsis et al.l ()2007l ) puts the limit at a 
more conservative 270 K, Kepler-69c would require an 
albedo greater than 0.5 to satisfy this requirement. On 
Fi gure [H we have plotted the habitable range suggested 
bv lKastiii^ (|2011f ). We have also plotted points corre- 
sponding to Venus, Earth and Mars. Kepler-69c likely 
falls within the ha bitable zone fo r moderate albedos if 
we use the limit of lKastina (|2011| ) although it would re- 
quire a fair ly high albedo to be habitable according to 
the limit of ISelsis et all (|2007[ ). 

The above analysis assumes that Kepler-69c has a cir- 
cular orbit. From our MCMC transit analysis we find 
Kepler-69c may be highly eccentric ~ the 3-cr upper limit 
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is e = 0.79. The apastron distance is 0.73 

which distance Kepler-69c would have an equilibrium 

temperature of 282^3^ K. However, Kcplcr-69c would 

have an equilibrium temperature of 322i2i K at perias- 
tron. If the planet is indeed on an eccentric orbit it may 
experience high seasonal temperature variations. How- 
ever, a Kepler-69c is most likely on a reasonably circular 
orbit and therefore will probably not go through dra- 
matic freczc-thaw cycles. 

7. CONCLUSIONS 




0.4 0.0 

Bond albedo 

Fig. 3. — The equilibrium temperature of Kepler-69c as a func- 
tion of Bond albedo. The blue shaded region shows the l-tr range 
assuming / = 1.0 where uncertainties were calculated from the 
MCMC chains. We expect / to be close to unity if Kepler-69c has 
an atmosphere similar to that of Earth. Venus, Earth and Mars 
are plot t ed for comparison. The habitable zone as suggested by 
IKastind I I2011I ) is shaded in red and bounded by the dashed lines. 
Kepler-69c falls near or within the habitable zone with reasonable 
assumptions of its albedo. 

In this paper we report the discovery of two transit-like 
signals in the light curve of Kepler-69 which we show are 
bona fide planets orbiting a Sun-like star. While we are 
confident at the 99.9% level that Kepler-69b orbits the 
target star, there is a slight possibility (< 0.9%) that 
Kepler-69c orbits a physical companion star to Kepler- 
69. 

With a l-cr uncertainty on the stellar radius and hence 
the radius of Kepler-69c of 17% we cannot conclusively 
say whe ther the planet is ro cky or has a high volatile 
content (jValencia et al.ll2007f ). 

The temperature of Kepler-69c falls close to its host 
star's habitable zone but a degree of inference must be 
made as to the planet's albedo in order to determine an 
equilibrium temperature. With reasonable assumptions 
it is possible that liquid water could exist in significant 
quantities on the surface of this planet. 

Kepler-69c represents the first discovery of a super- 
Earth in the habitable zone of a Sun-like star. The work 
represents a progressive step on the road to detecting the 
first truly Earth-like planet orbiting a star like our Sun. 

This paper includes data collected by the Kepler mis- 
sion. Funding for the Kepler mission is provided by the 
NASA Science Mission Directorate. Some Kepler data 
presented in this paper were obtained from the Mikulski 
Archive for Space Telescopes (MAST) at the Space Tele- 
scope Science Institute (STScI). STScI is operated by 
the Association of Universities for Research in Astron- 
omy, Inc., under NASA contract NAS5-26555. Support 
for MAST for non-HST data is provided by the NASA 
Office of Space Science via grant NNX09AF08G and by 
other grants and contracts. This research has made use 
of the NASA Exoplanet Archive, which is operated by 
the California Institute of Technology, under contract 
with the National Aeronautics and Space Administra- 
tion under the Exoplanet Exploration Program. We used 
data from the XJK IRT Wide Field Camera (WFCAM, 
iCasali et al.ll2007l ) and a photometric system described 
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in lHewett et al.l (I2006D . The pipeline proces s ing an d sci- 
ence archive are described in iHamblv et al.l (|200^. DH 
is supported by appointment to the NASA Postdoctoral 



Program at Ames Research Center, administered by Oak 
Ridge Associated Universities through a contract with 
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